Vibration-induced jamming transition in granular media 
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The quasi-static frequency response of a granular medium is measured by a forced torsion oscillator 
method, with forcing frequency f p in the range 10 -4 Hz to 5 Hz, while weak vibrations at high- 
frequency f s , in the range 50 Hz to 200 Hz, are generated by an external shaker. The intensity of 
vibration, T, is below the fluidization limit. A loss factor peak is observed in the oscillator response 
as a function of V or f p . In a plot of In f p against 1/F, the position of the peak follows an Arrhenius- 
like behaviour over four orders of magnitude in f p . The data can be described as a stochastic hopping 
process involving a probability factor exp(— Tj/T) with Vj a / 3 -dependent characteristic vibration 
intensity. A / s -independent description is given by exp(— Tj/r), with Tj an intrinsic characteristic 
time, and r = T n /2nf s , n=0.5-0.6, an empirical control parameter with unit of time, r is seen as the 
effective average time during which the perturbed grains can undergo structural rearrangement. The 
loss factor peak appears as a crossover in the dynamic behaviour of the vibrated granular system, 
which, at the time-scale 1/ f p , is solid-like at low F, and the oscillator is jammed into the granular 
material, and is fluid-like at high T, where the oscillator can slide viscously. 
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I. INTRODUCTION 

In a granular medium at rest the grains can be disposed 
in an enormous number of different configurations. A 
weak external disturbance, but powerful enough to over- 
come locally the friction force between two grains, allows 
the granular systems to rearrange and to switch between 
these "blocked" configurations. The macroscopic behav- 
ior of a weakly disturbed granular medium is, therefore, 
essentially controlled by statistical properties of such 
transitions. It is of great interest to study this kind of 
problem, since it may be a prototype of slow dynamics 
behaviors observed in other physical systems Q. The 
slow dynamics of weakly disturbed granular media has 
been evidenced by the classical compaction experiments 
of Knight et al. \ 3 Another approach is the experiment 
of Albert et al. ]3], in which a large solid object is pulled 
slowly through a granular medium. The motion of the 
object appears as resisted by chains of jammed parti- 
cles Q 0, which support compressive stress. Beyond an 
elastic regime at very small pulling force, the macroscopic 
motion of the object is a succession of stick-slip events, 
where compressive stress is continuously built up in par- 
ticle chains, and abruptly released. At these successive 
unjamming events the system switches between blocked 
configurations. 

In this paper we implement a novel experimental 
method to study the problem: we exploit a forced tor- 
sion oscillator j^] immersed in the granular medium, and 
in presence of external weak vibration, as shown in Fig. 
1. We use a dynamic method, which provides more infor- 
mation than the simple increase toward the unjamming 
threshold. In fact, the amplitude of the angular displace- 
ment of the oscillator increases sharply when the unjam- 
ming threshold forcing torque amplitude is approached, 



and at the same time the angular displacement lags be- 
hind the sinusoidal torque. This "phase lag" is deter- 
mined by an energy dissipation which occurs in the gran- 
ular medium when grains start to slip one against the 
other. The dynamic method gives access to both elas- 
tic and dissipation parameters of the granular material 
during the slow dynamics. 



II. EXPERIMENTAL 

In the experiment, we hold a granular medium at a 
given high-frequency vibration intensity, quantified by 
the normalized acceleration T = a s ajg/g, with a s and 
f s = lu s /2tt the amplitude and frequency of the verti- 
cal sinusoidal vibration, g the acceleration of gravity. At 
the same time, we measure the complex frequency re- 
sponse, G, of the granular medium (or the susceptibility 
X = G~ x ) by a low- frequency forced torsion oscillator M, 
at the forcing frequency f p — w p /27r, with f p <C ,f s . In 
the oscillator method (see Fig. 1), the rotating probe 
of the oscillator is immersed at a depth L into a large 
metallic bucket (height 96 mm, diameter 94 mm) filled 
with glass beads of diameter d = 1.1 ± 0.05 mm with 
smoothly polished surfaces. The probe is covered by a 
layer of beads, glued on by an epoxy, and its effective 
radius is R e . All data presented here are obtained with 
L = 20 mm and R e ~ 2 mm. We perform dynamic ex- 
periments: the oscillator is forced into torsion oscillation 
by a torque T (t) — T Q exp(iu)pt) of frequency in the range 
10~ 4 Hz to 5 Hz, and the angular displacement, 0(t), is 
optically detected. An analyzer measures the complex 
frequency response of the oscillator, given by G = T/6. 
Typically we record the argument, arg(Gi), and the abso- 
lute modulus, |Gi|, of the first harmonic, as a function of 
either T , T, or f p . We report the quantity tan[arg(Gi)], 
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which for a linear system coincides with the loss factor. 
The oscillator, when not immersed, can be assumed elas- 
tic, with T = G p 9, where G p = 18 x 1CT 3 N m/rad is the 
torsion constant of the suspension wires. Notice that in 
this work the maximum displacement of a point at the 
surface of the 2 mm probe is of the order of 0.1 mm, i.e., 
much smaller than the glass bead diameter. 

An accelerometer provides a precise measurement of 
r, which can be varied from 2 x 10~ 3 to above 1. The 
minimum value of T is limited by the accelerometer sen- 
sitivity. We vary T by changing a s at fixed f s , while f s 
is selected in the range 50 Hz to 200 Hz. The whole sys- 
tem is placed on an anti-seismic table. Moisture-induced 
ageing effects |J, and interstitial gas effects ||, are 
not observed for the large bead size used here, and mea- 
surements are performed at uncontrolled ambient air. In 
order to control compaction effects [|| , all measurements 
are taken in the same conditions, e.g., starting from a 
granular material shaken at high T and low f s for several 
minutes. Compaction effects are apparently negligible in 
the time scale of the experiments for f p > 0.01 Hz, but 
may be present in the data at very-low frequency. 

III. RESULTS 

A typical experimental result is reported in Fig. 2, 
which shows tan[arg(Gi)] and \G\\, measured as a func- 
tion of the amplitude of the applied torque To, for differ- 
ent T. With the vibrator off, that is for T < 2 x 10~ 3 , 
the response is similar to the one reported previously Q , 
with a typical loss factor peak at a torque denoted TJ 
and a modulus step between two levels denoted Gj am and 
G p . The dependence of the loss peak on the geometrical 
parameters of the experiment is summarized [0 by the 
empirical relation T$ ex fi s L 2 R 2 , with /i s the coefficient 
of static friction between the glass beads. In this "zero 
temperature-like" conditions, the loss factor peak can be 
easily explained: at very low applied torque, T <C T *, 
the oscillator probe is jammed into the granular material, 
and only elastic deformations arise, resulting in a purely 
elastic dynamic response, with a negligible loss factor 
and a constant absolute modulus Gj am - By increasing 
the torque amplitude, the oscillator probe unjams as the 
local force between a pair of glass beads somewhere in the 
medium becomes large enough for the two beads to slip 
one against the other, dissipating energy by solid friction. 
The maximum ratio of dissipated over furnished energy, 
that is a maximum of the loss factor, arises at T *, which 
can be seen as the average torque at which the oscillator 
probe unjams. At high torque amplitude, To 3> Tg , the 
oscillator slides almost freely into the granular medium, 
and the absolute modulus tends to the torsion constant 
of the suspension wires, G p . 

With the vibrator on, one expects that the external 
vibration facilitates the unjamming of the oscillator. By 



increasing T, the modulus |Gi| decreases monotonically, 
while tan[arg(Gi)] first increases and then decreases, go- 
ing through a maximum, as clearly visible in Fig. 2. 
The fact that the external vibration drives the system 
through a maximum in the loss factor is evidence that 
the vibration-induced fluctuations can unjam the oscilla- 
tor probe. Moreover, below T * the response is essentially 
independent of To , i.e., there is a linear regime. The lin- 
earity is confirmed also by a negligible high harmonics 
signal (not shown) for all T. 

The general behavior in the linear regime is better ren- 
dered by Fig. 3 which shows the previous maximum 
in the loss factor as characteristic "jamming" peaks ob- 
served as a function of T for various forcing frequencies 
f p . In Fig. 3 we obtain that for the low-torque amplitude, 
selected in the linear regime, i.e., To T *, and at low T, 
the applied torque alone is unable to unjam the oscillator 
probe and the response is elastic. However, by increasing 
r, unjamming is induced by the external vibration, and 
the response displays a loss factor peak. The data shown 
in Fig. 3 are collected by decreasing T, but no difference 
is observed in following runs if T is successively increased, 
decreased and so on, as shown in Fig. 4 for one of the 
curves of Fig. 3. We say that the response is " reversible" , 
although at a mesoscopic level energy is continuously dis- 
sipated. The loss factor peak can be seen as the crossover 
between two different behaviours in the dynamics of the 
vibrated granular system: at the time-scale set on by the 
forced oscillator, i.e., 1/ f p , the granular system appears 
solid-like at low-r, while it appears fluid-like at high-r. 
It is a kind of glass, or jamming transition pp[ | , where 
the oscillator gets stuck in the glassy granular medium. 

Of course, since the "jamming" peak in Fig. 3 shifts 
with f s , the same peak can be observed as a function of 
the forcing frequency f p , as shown in Fig. 5, for vari- 
ous r. At high forcing frequency (but still f p <C / s ), the 
applied torque alone is unable to unjam the oscillator 
probe and the response is elastic, with negligible loss fac- 
tor and modulus Gj am . However, by decreasing f p , i.e., 
by increasing the time scale of the probing oscillator, the 
response evolves toward the one of the unjammed oscil- 
lator, with modulus G p . 

From the shift of the previous "jamming" peaks with 
fp or r, an Arrhenius-like semilogarithmic plot can be 
obtained, as shown in Fig. 6. The data points, for a 
given f s , obey an exponential behavior over four decades 
in frequency. This is strong evidence for the underly- 
ing unjamming process to be a statistical, activated-like 
hopping process and that some of the usual statistical 
concepts of thermal systems can be extended to a vi- 
brated granular material. A first, obvious approach con- 
sists of formally writing the rate of the hopping process as 
R = cxp(— Tj/T), with Tj a characteristic normalized 
acceleration at which unjamming occurs, vq an attempt 
frequency, and T the vibration intensity, playing the role 
of a temperature-like parameter. In the linear regime, 
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a peak in the loss factor is expected to arise when the 
forcing frequency matches the hopping rate, i.e., when 
ojp = R, and Tj appears as the "slope" of a straight line 
in a plot of \nR against 1/T; however, then Tj depends 
on f s (see Fig. 6), which means that Tj is not an intrinsic 
property of the granular material. 

To overcome this difficulty, we search a "scaling" of 
T and f s which eliminates the f s dependence. We find 
that as a function of the inverse of T n /ui s , with n = 0.5, 
the data for different f s have almost the same " slope" , as 
shown in Fig. 7. Hence, we write the hopping probability 
per unit time as R = vq exp(— tj/t), with tj sl character- 
istic time of the unjamming process, and r = r 1 / 2 juj s a 
parameter which has the unit of time |j| . (Fitting values 
are given in Fig. 7) Alternatively (Fig. 8), we find also 
that a "scaling" of the form r = T n /oj s , with n — 0.57, al- 
most collapses the data of Fig. 6 over the same straight 
line. This empirical definition is appealing since both 
the parameters r and Tj we introduce, are independent 
of the vibration frequency f s . The unjamming time Tj 
experimentally is an intrinsic parameter of the granular 
system, which for a given L and R e is likely to depend 
on "mesoscopic" parameters such as the grain size and 
shape, and on "microscopic" parameters controlling the 
nature of the contact forces between grains. For n = 0.5 
the externally controlled parameter r appears as a typical 
time of the vibrated granular system in the gravitational 
field: V 1 / 2 /uj s — (a^/g) 1 / 2 is the time of flight of a body, 
initially at rest, falling for a distance a s . (It is also the 
period of a simple oscillator with thread a s and freely 
swinging.) Notice that r is not a kinetic energy- type 
temperature, as defined for vigorously vibrated gas-like 
granular phases; indeed, for weakly excited granular sys- 
tems, configurational statistics on slow degrees of free- 
dom can be decoupled from fast kinetic aspects fl |l(]] 

IV. DISCUSSION 

That the rate of our " activated- like" process is better 
given by a probability factor involving the ratio of char- 
acteristic times, exp(— Tj/r), and not, e.g., by the ratio of 
vibration intensities, exp(— Tj/T), is not surprising since 
the fundamental phenomenon is the rate of energy dis- 
sipation. What is the mesoscopic nature of the hopping 
processes at the length scale of the glass beads? Since we 
observe an elastic regime at low T, we conclude that the 
oscillator probe is completely jammed and no dissipative 
events (i.e., no slipping events) arise. We can suppose the 
system to oscillate elastically around one unique blocked 
configuration. In the picture, below the loss factor peak, 
the external high-frequency vibration propagates in the 
system as elastic fluctuations only. Such elastic fluctu- 
ations are non-dissipative vibrations of the force-chain 
network which holds the oscillator probe in place. By in- 



creasing r, a large number of increasingly energetic elas- 
tic fluctuations arise (or increasingly longer if we focus on 
the parameter r = T n /uj s ), and the force-chain network 
can be seen as exploring different elastic configurations 
until a critical configuration is reached. A critical config- 
uration is such that, somewhere, the local friction force 
between two glass beads is overcome and slipping arises, 
momentarily unjamming the oscillator probe. The sys- 
tem can switch to another blocked configuration. 

A single slipping event possibly triggers a large-scale 
non-elastic rearrangement of the beads, that is an inter- 
nal micro-avalanche involving two or more beads. (We 
emphasize that a massive object moving into a granular 
medium can introduce local fluidization. The inertia ac- 
quired by the object during a slip may be large enough 
to overcome the resisting force of the granular medium, 
and the object moves further by successive failure, or "in- 
ertial" fluidisation, of the resisting grains arrangement. 
This effect can be reduced by immersing the object deep 
enough, so that the resisting force is larger than the iner- 
tial force.) Afterwards, the oscillator probe gets jammed 
by a new force-chain network in a slightly different posi- 
tion. The macroscopic slow rotation of the oscillator is a 
sequence of stick-slips, where large fluctuations causing a 
slip are rare events if compared to the numerous elastic 
fluctuations. Hence, the dynamics is controlled by the 
extreme fluctuations in the force-chain network, even if 
elastic, or under-critical fluctuations occur in much larger 
number. Since a slipping event involves inelastic micro- 
scopic processes at the interface between grains, such 
as plastic and viscoplastic deformation, fatigue, surface 
fracture, blow out of capillary bridges, and other forms 
of localized dissipative processes |L5j , a slipping event re- 
quires a minimum finite time to occur. Let this time be 
Tj. The "thermal" time r can be seen as the average 
time- window during which the grains have some freedom 
to rearrange their position and, possibly, reach a criti- 
cal slipping configuration and unjam. As a consequence, 
unjamming is determined by the occurrence probability 
of a window time T a greater as, or equal to Tj. Even 
though we do not know the probability distribution of 
t q , according to extreme order statistics theory (T(| 
Pq| pl| , we can speculate the occurrence probability of 
unjamming events to be exp(— Tj/r). This gives the rate 
of the extreme fluctuations, namely R = vq exp(— Tj/r). 

In summarizing, we observe a peak in the loss factor 
as a function of the empirical control parameters r or T. 
The peak can be viewed as a crossover, at the time-scale 
l// p set on by the forced oscillator, in the dynamics of 
the vibrated granular medium: such a crossover sepa- 
rates a "low-temperature" (short r, or small T) solid-like 
behaviour where the oscillator probe is jammed in the 
granular medium, from a " high-temperature" (long r, or 
large T) fluid-like dynamic behaviour. This crossover fol- 
lows an Arrhenius-like form in the In/ vs. 1/T plane, 
reminiscent of the mechanical response of usual glass- 
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framing materials. 
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Figure captions 

FIG. 1. Sketch of the forced torsion oscillator im- 
mersed at a depth L into a granular medium of glass 
beads. A single layer of glass beads is glued to the os- 
cillator probe, and the effective radius is R e . The con- 
tainer, filled with the granular material, is shaken verti- 
cally by a vibrator at the intensity of vibration, T, below 
the fiuidization limit. The method provides a measure of 
the complex frequency response of the granular medium 
while the vibrator mimics " thermal" fluctuations, ^sus- 
pension wires; 2=permanent magnet; 3=external coils; 



4=mirror; 5=probe; 6=vibrator. 

FIG. 2. Oscillator frequency response, at f p = 1 Hz, 
as a function of the amplitude of the applied torque To, 
for different vibration intensities T, for f s = 200 Hz. (a) 
The loss factor, given by tan[arg(Gi)], versus To. The 
position of the peak obtained at T < 2 x 10~ 3 (i.e., with 
the vibration off) is denoted Tg. (b) The modulus |Gi| 
versus T). The two extreme levels on the curve obtained 
at T < 2 x 10 -5 are denoted G p and Gj am respectively. 
For T) <§C T *, the response is independent of To , i.e., 
there is a linear regime, confirmed also by a negligible 
high harmonics signal (not shown) for all T. 

FIG. 3. Oscillator frequency response as a function 
of the vibration intensity T, with f s = 200 Hz, for dif- 
ferent forcing frequencies f p of the oscillator, at a given 
low torque amplitude To = 3.2 x 10~ 5 N m selected in 
the linear regime, i.e., To -C T *. (a) The loss factor 
tan[arg(Gi)], versus T. (b) The modulus |Gi| versus T. 
For each f p a peak with a maximum at a vibration in- 
tensity r* can be seen. The data shown in Fig. 3 are 
collected by decreasing T, but no difference is observed 
in following runs if T is successively increased, decreased 
and so on (see Fig. 4.). 

FIG. 4. Similar to Fig. 3, for f p = 1 Hz, but for T 
successively decreased, increased, decreased and so on. 

FIG. 5. Oscillator frequency response as a function of 
the forcing frequency f p , for different vibration intensity 
r, with f s = 200 Hz, at a given torque amplitude To — 
3.2 x 10 -5 N m. (a) The loss factor tan[arg(Gi)], versus 
f p . (b) The modulus |Gi| versus f p . One can see for 
each r a peak with a maximum at a frequency f*. The 
data shown are collected by decreasing f p . A Debye peak 
of equation Coj p t c /(1 + uJ p r^) with G = 1.9 and r c = 
= 0.6 is shown in (a). As a function of the frequency, 
the shape of the Debye peak is independent of the exact 
definition of the temperature-like parameter entering the 
rate R. The observed "jamming" peaks are much larger 
than the pure Debye peak, suggesting that the underlying 
dynamics is glassy in nature. 

FIG. 6. The semilogarithmic Arrhenius-like plot re- 
porting the forcing frequency f p versus 1/r of "jamming" 
peaks similar to the ones in Figs. 3 and 5 (filled symbols 
from measurements vs. T; open symbols from measure- 
ments vs. f p ), for different vibration frequencies f s . The 
data are fitted (dashed lines) by 27r/ p = v exp(— Tj/T), 
which gives v «70 Hz and Tj «0.014 for f s = 50 Hz, 
u «66 Hz and Tj «0.042 for f s =100 Hz, and v «49 Hz 
and T 3 wO.14 for f s =200 Hz. 

FIG. 7. The forcing frequency f p versus the inverse of 
the empirical control parameter 1/t, i.e., versus oj s /T n , 
with n — 1/2. The data are fitted (plain and dashed 
lines) by 27r/ p = v> exp(— Tj/r), which gives vq wl336 Hz 
and Tj w 1.3 x 10~ 3 s for f s = 50 Hz, v «1465 Hz and 
tj « 1.2 x 10~ 3 s for f s = 100 Hz, u w2698 Hz and 
Tj w 1.2 x 10~ 3 s for f s = 200 Hz. The straight lines have 
almost the same "slope" Tj. The average is (Tj) — 1.26 x 
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1CP 3 s. vq is seen as a natural vibration frequency of the 
granular medium. Considering the present precision of 
the data, no clear relationship between v and f s can be 
found, even though vq increases as f s increases. 

FIG. 8. Similar to Fig. 7, but with n = 0.57. 
The data are fitted (plain and dashed lines) by 2nf p — 
i/o exp(— tj/t), which gives v$ «648 Hz and tj « 8.0 x 
10~ 4 s for f s = 50 Hz, v «683 Hz and Tj « 7.8 x 10~ 4 s 
for f s = 100 Hz, v wl019 Hz and Tj « 8.7 x 10~ 4 s for 
f s = 200 Hz. 
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